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ABSTRACT 

The formation of dark matter halos tends to occur anisotropically along the filaments 
of the Cosmic Web, which induces both ellipticity-eUipticity (EE) correlations between 
the shapes of halos, as well as ellipticity-direction (ED) cross-correlations between 
halo shapes and the directions to neighboring halos. We analyze the halo catalogue 
and the semi-analytic galaxy catalogue of the recent Millennium Run Simulation to 
measure the EE and ED correlations numerically at four different rcdshifts (z — 0, 
0.5, 1 and 2). For the EE correlations, we find that (i) the major-axis correlation is 
strongest while the intermediate-axis correlation is weakest; (ii) the signal is significant 
at distances out to 10/i^^Mpc; (iii) the signal decreases as z decreases; (iv) and its 
behavior depends strongly on the halo mass scale, with larger masses showing stronger 
correlations at large distances. For the ED correlations, we find that (i) the correlations 
are much stronger than the EE correlations, and are significant even out to distances 
of 50/i~^Mpc; (ii) the signal also decreases as z decreases; (iii) and it increases with 
halo mass at all distances. We also provide empirical fitting functions for the EE and 
ED correlations. The EE correlations are found to scale linearly with the linear density 
correlation function, ^(r). While the ED cross-correlation is found to scale as S,^^^{r) 
at large distances beyond 10/i~^Mpc. The best-fit values of the fitting parameters for 
the EE and the ED correlations are all determined through x^-statistics. Our results 
may be useful for quantifying the filamentary distribution of dark matter halos over 
a wide range of scales. 

Key words: methods:statistical - cosmology:theory - galaxies: clustering - galax- 
ies:halos - large-scale structure of Universe 



1 INTRODUCTION 

One of the most striking features of the Universe is that 
the observed distribution of galaxies on large scales shows a 
web-like filamentary pattern, which is often called the "Cos- 
mic Web" . Recent large N-body simulations of the cold dark 
matter cosmology demonstrated vividly the geometric rich- 
ness of the filamentary web that spatially connects the dark 
matter halos, and which directly relates to the structure 
seen in the galaxy distribution. One of the most fundamen- 
tal tasks in cosmology is thus to establish a physical model 
for the filamentary cosmic web and to quantitatively explain 
its global properties. 

The existence of the filamentary web was originally 
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predicted by the to p-down scenario of the hot dark mat- 
ter (HDM) model (IZel'dovichl ligTOD . If cosmic structures 
form through top-down fragmentation, then one- and two- 
dimensional collapse of matter would naturally lead to the 
formation of sheet-like and filamentary structures on large 
scales. Therefore, it was regarded first as a mystery why and 
how the filamentary web came into being also in a cold dark 
matter (CDM) dominated universe. 

A breakthrough was made by iBond et al.l ||1996D who 
developed a cosmic web theory that can explain the pres- 
ence of a filamentary web in the CDM cosmogony. This the- 
ory explains that the filamentary web can occur naturally 
in the CDM dominated universe due to the coherent na- 
ture of the primordial tidal field. The filamentary web is in 
fact a manifestation of the primordial tidal field sharpened 
by nonlinear effects. The cosmic web theory has provided a 
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standard framework within which the formation of cosmic 
large-scale structure can be qualitatively understood. Yet, 
it is still quite difficult to describe the cosmic web quantita- 
tively both in theoretical and in observational terms. The- 
oretically, the inherent anisotropic nature and geometrical 
complexity of the cosmic web makes it complicated to fully 
characterize its statistical properties. Observationally, it is 
hard to trace the filamentary structures from observational 
data, since there is no well-established way to identify them. 

In spite of these difficulties, various methodologies and 
algorithms have already been suggested to quantify the fila- 
mentary structures: Higher-order N-point statistics has been 
used to desc ribe the anisotropic matter distribution in a 
cosmic web l|Croton et al.ll2004l : [Kulkarni et al.l 120071 '): the 
percolation statistics was used to charac terize the filamen- 
tary shapes of the large-scale structures l|Sahni et al.lll997l : 
[Shandarin & Ycss 1998); the skeleton formalism has been 
developed to extract the filam entary structures fr om a three 
dimensional density field (e.g. JSousbie et al]|2007l . and refer- 
ences therein); the Minimal-Spanning- Three algorithm has 
been introduced to find the basic structural elements of the 
cosmic web ( Colborg 2 007i) . 

Although the above methods are quite useful for de- 
termining the overall filamentary structure of the cosmic 
web, these approaches are largely phenomenological with- 
out accounting for the physical mechanism for the forma- 
tio n of the cosmic w eb. According to the theory proposed 
bv lBond et al.l (|l996h . the filamentary web originated from 
the large-scale coherence of the primordial tidal field and 
its sharpening by nonlinear effects during structure growth. 
Part of this nonlinear sharpening effect arises from the grav- 
itationally driven merging of halos and the infall of matter, 
which preferentially occurs along the most prominent fila- 
ments. This increases the anisotropy in the halo clustering 
and thus sharpens the filamentary web. 

To describe the cosmic web quantitatively in terms of its 
underlying physical principles it will be necessary to account 
for the effects of the tidal field and the anisotropic merging 
along filaments. The tidal field causes intrinsic alignments 
of the principal axes o f the dark halo s in the cos mic web 
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anisotropic merging induces elongati on of the ma.jor axes of 
the h alos along prominent filaments (jWestll 19891 : IWest et al.l 
I1991] ). As a result, there exist spatial correlations between 
the halo ellipticities (EE correlations), and cross-correlations 
between the halo ellipticities and the large-scale density field 
(ED cross-correlations). Hence, one can view the observed 
filamentary web as a large-scale manifestation of the EE and 
ED correlations, which are in turn induced by the effects of 
the tidal field and the anisotropic merging. 

The goal of this paper is to quantify the filamentarity 
and the typical scales of the cosmic web in terms of the EE 
and ED correlations. This is also highly important for assess- 
ing to what degree these correlations can systematically bias 
weak gravitational lensing mass reconstructions and cosmo- 
logical parameter estimates based on cosmic shear measure- 
ments. In fact, the ED cross-correlations have become a hot 
issue in the weak lensing community, since it has been real- 
ized that t hey could mimic weak lensing signals at a signif- 
icant level (iHirata fc Seliakll2004 iMandelbaum et all [200^ : 
iHirata et al.ll2007l '). 



Several authors have already studie d num erically the 
EE and ED correlations. iHopkins et al.l (|2005l ) examined 
the evolution of cluster a lignments by using a large-scale N- 
body simulation. lAltav et al. (2006) have shown from high- 
resolution N-body data that the alignments of clusters are 
strongly related to the existence of the connecting filaments. 
iHevmans et al.l (|2006l ) explored the possible correlations be- 
tween the weak lensing shear and the orientation of fore- 
ground galaxies by analyzing N-body simulation. Here in 
our work, we measure both the EE and ED correlations at 
large distances by using the Millennium dataset and describe 
quantitatively their scalings with distance, mass, as well as 
redshift. 

The organization of this paper is as follows. In Section 2, 
we describe the N-body dataset we use and explain how we 
measure halo ellipticities from the N-body simulation and 
its associated galaxy catalogue. In Sections 3 and 4, we re- 
port numerical detections of the EE correlations and the ED 
cross-correlations, and examine how the signals depend on 
distance scale, redshift and halo mass. We also provide useful 
fitting formula for them. Finally, in Section 5, we summarize 
the results and discuss the implications of our work. 



2 SIMULATION DATA AND METHODOLOGY 

Our analysis is based on the halo catalog and the semi- 
analytic galaxy catalog from the recent high-resolution Mil- 
lennium Simulatio'^ which followed 10^" da rk matter par- 
ticles in a ACDM concordance cosmology (jSpringel et al.l 
2005). The size of the periodic simulation box is 500 /i~^Mpc 
and each dark matter particle in the simulation has a mass 
of 8.6 X lcPh~^ Mq. The basic cosmological parameters of 
the simulation were chosen as Q,m = 0.25 (the mass den- 
sity); = 0.75 (the vacuum energy density); h = 0.73 
(the dimensionless Hubble constant); erg = 0.9 (the linear 
power spectrum amplitude); and = 1 (the slope of the 
primordial power spectrum). 

As part of the analysis of the Millennium run, halos of 
dark matter particles were first identified with the standard 
friends-of-friends (EOF) algorithm, and then decomposed 
into gravi tationally bound sub halos using the SUBFIND al- 
gorithm (jSpringel et al.l [20011 ). Based on detailed merger 
trees constructed for the subhalos, the halos were then popu- 
lated with luminous galaxy models using se mi- analytic simu - 
lations of the galaxy formation process (Cr oton et aLlbOOd ) . 

We here use the spatial distribution of subhalos and 
galaxies to characterize the shape of FOF halos, in analogy 
to th e procedure ap plied to observational galaxy surveys 
(e.g., iMei et al.1 [20071 '). For each FOF halo, we locate the 
satellite galaxies belonging to it. Then, we measure their 
tensor {hj) of second order mass moments as 



(1) 



where ma is the luminosity (or, equivalently the stellar 
mass) of the a-th galaxy and Xa is the position of the a-th 
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galaxy measured from the center of the mass of the satel- 
lite galaxies. We restrict our analysis to FOF halos massive 
enough to contain more than five substructures. By diago- 
nalizing lij, we determine the three principal axes (major, 
intermediate, and minor axes) of lij . This allows us to mea- 
sure the correlations between the three axes of the FOF 
halos as a function of separation. The results of our mea- 
surements are presented in detail in Section 3. 

Before turning to our results, it is worth to discuss how 
our methodology relates to other, previously applied meth- 
ods to characterize the shape of halos. Note that we here do 
not use all dark matter particles of a FOF halo to measure 
lij. Instead, we only use the satellite galaxies (or substruc- 
tures) as tracers of the shape. In general, measuring the 
shape of a halo is a somewhat ambiguous issue, where a 
number of different strategies have been applied in the liter- 
ature, but no generally acce pted sta ndard procedu re exists 
(see e .g. the discussion in Spririgel et al. 2004; AUg ood et al.l 
I2OO6I ). Part of the ambiguity in measuring halo shape stems 
from the fact that one cannot delineate the outer boundary 
of a halo in a clear-cut way. If all particles belonging to a 
FOF halo are used to measure lij, then the ellipticity of a 
halo may be overestimated because of the large weight of 
the most distant points on the major axis, while in contrast, 
if only those particles within a certain spherical radius are 
used to measure hj , then the ellipticity of the halo is likely 
underestimated. 

We use satellite galaxies (or substructures) to measure 
lij , because this approach mimics observationally accessible 
procedures. We expect that this definition should give halo 
shapes similar to those measured with all dark matter par- 
ticles, as substructure and gal axy density are tracers of th e 
dark matter distribution (e.g.. lAgustsson fc BrainerdllioO^ . 
We will explicitly test this below. Note that we focus on 
quantifying the filamentary web induced by the anisotropic 
orientations of the halo ellipticities; we are not really inter- 
ested in the magnitude of the shape distortion itself. This 
means that we are less sensitive to the details of measur- 
ing halo shape compared with attempts to quantify the axis 
ratio of the halo shape. 

To examine to what extent different measuring methods 
yield different halo ellipticities, we carry out a simple test: 
Using a total of 227 FOF halos from the 'milli'-Millennium 
simulation, which is smaller test run of the m ain Millennium 
simu lation with a box size of 62.5 /i~^Mpc jSpringel et al.l 
I2OO5I ). we measured the halo shapes with three different 
methods: (A) using galaxies weighted by luminosity; (B) us- 
ing dark matter substructures weighted by mass; and (C) 
using all dark matter particles. Then, we calculate the dis- 
tribution of the angles 6 between the principal axes of the 
halo shapes determined by these three methods. Figure [T] 
compares methods A and C (top), and methods B and C 
(bottom), in both cases plotting the histogram of the an- 
gles between the halos' major, intermediate, and minor axes 
(left, middle, and right panels, respectively). As can be seen, 
there is a strong peak ai — 0, demonstrating that the halo 
principal axes obtained by the three different methods A, B 
and C are strongly correlated with one another. This corre- 
lation is particularly robust for the major axis, which defines 
the primary orientation of the predominantly prolate halos. 
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Figure 1. Distributions of the angles between the major, inter- 
mediate, and minor axes of the halos (left, middle, and right, re- 
spectively) from the milli-Millennium simulation, as determined 
by three different methods A, B, and C, which are based on the 
satellite galaxies, subhalos, and all particles belonging to the ha- 
los, respectively. The top-panels show the results from a compar- 
ison between the methods A and C, while the bottom panels give 
a comparison between the methods B and C. 



3 THE HALO ELLIPTICITY-ELLIPTICITY 
CORRELATION 

3.1 Definition 

We define the EE correlation function, ri{r), as 

77(r)^(|e(x).e(X + r)|^)-i (2) 

where e = (e^) represents the normalized eigenvector of a 
halo with unit magnitude. In eq.((2} the constant 1/3 is sub- 
tracted since the first average term will yield 1/3 in case that 
there is no EE correlation. Now that the three eigenvectors 
of each halo in the Millennium catalogs are all determined 
by the method described in Section 2, one can measure the 
EE correlations of the major (771), intermediate (7711) and 
minor (77111) principal axes separately as a function of the 
comoving distance r. 

Basically, for each halo pair at a given redshift, we mea- 
sure the separation distance r between the halo centers and 
calculate the squares of the dot products of the normalized 
eigenvectors of two halos. Then, we bin the radial distance 
r and calculate the mean values of |e(x) -ejx-l-r)!^ averaged 
over those halo pairs whose separation distances belong to 
a given bin, subtracting 1/3 from it. We perform this pro- 
cedure at 2 = 2, 1, 0.5 and 0. 



3.2 Evolution with redshift 

Figure [5] plots the numerical results on the EE correlations 
measured at z = 0, 0.5, 1 and 2 in the top-left, top-right, 
bottom-left and bottom-right panels, respectively. In each 
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Figure 2. EE correlations of the halo major, intermediate, and 
minor axes (solid, dashed, and long dashed lines) at four different 
redshifts: 2 = 0, 0.5, 1, and 2 (top-left, top-right, bottom-left, and 
bottom-right panel, respectively). 

panel, the solid, dashed and long dashed lines represent 
rjiir), 7)11 (r) and r)iii{r), respectively. The dotted line cor- 
responds to the case of no correlation. As can be seen, the 
major-axis correlations are strongest and the intermediate- 
axis correlations are almost zero at all redshifts. 

To see the behaviors of the EE correlations at large dis- 
tances, we plot 771 (r) and r)iii{r) as solid dots with errors 
CTj, in the logarithmic scale in Figs [3] and |4j respectively. In 
each panel the solid line represents the fitting model (see 
Section 3.4). For the estimation of (Jr,, we divide the simula- 
tion volume into eight subvolumes each of which has a linear 
size of 250/i~^Mpc and measure the EE correlations in each 
subvolume separately. The errors, cr,,, are calculated as the 
standard deviation between realizations. This estimation of 
errors accounts for both the cosmic variance and the Poisson 
noise. It is also found that there exist non-negligible correla- 
tions between radial bins at distance larger than 5/i~^Mpc 
(see Section 3.4). 

As can be seen, the EE correlations are strongest at 
2 = 2, and exist out to distances of 10/i~^Mpc. As 2 de- 
creases, the correlations tend to decrease monotonically at 
all distance scales, indicating that the directions of the halo 
major and the minor axes tend to become randomized as z 
decreases. This result is consisten t with the previous results 
obtained bv lHopkins et all (|2005l ). The numerical results on 
the EE correlations measured at 2 = are listed in Table [T] 



3.3 Variation with mass 

To study how the EE correlations scale with halo mass, we 
measure r?i(r) from two different mass bins with the mass 
threshold Mc = 1.42 x IO^^/i'^Mq at z = 0. Table [2] lists 
the numerical results and Figure [5] plots 771 (r) at z = mea- 



Figure 3. The EE correlations of the halo major axes at 2 = 0, 
0.5, 1 and 2. The errors, cr^, represent the standard deviation 
between 8 realizations that are obtained from the subdivision of 
the simulation box. The errors do not include the correlations be- 
tween radial bins. In each panel the solid line represents the fitting 
model (eqlSjl proportional to the linear density correlation func- 
tion, i,(r). In the bottom-right panel, the dashed line represents 
another fitting model proportional to ^^(r). 



Table 1. Numerical results on the EE correlations of the halo 
major (rji) and the minor axes (rjni) in logarithmic scale measured 
at z = 0. 



log[r/(/i-iMpc)] 


rii{r) X 10^ 


?7lll(r) X 10^ 


0.15 


1.00 ±0.69 


1.04 ±0.97 


0.45 


0.94 ±0.22 


0.89 ±0.38 


0.75 


0.58 ± 0.12 


0.44 ±0.12 


1.05 


0.22 ± 0.08 


0.10 ±0.06 


1.35 


0.09 ± 0.03 


0.03 ±0.02 


1.65 


0.02 ± 0.02 


0.01 ±0.02 



sured from the low-mass (M < Mc) and the high-mass bin 
(M > Mc) as solid dots in the top and the bottom panel, 
respectively. As can be seen, the EE correlations of the high- 
mass halos are stronger at all distances than that of the low- 
mass halos, which implies that the EE correlations increase 
as the halo mass increases. This finding is also consistent 
wi th the previous results obtained for the cluster alignments 
bv lHopkins et al.1 (|200d ). 

We have also measured the EE cross-correlations, r]c{r ), 
between the low-mass and the high-mass halos. Basically, we 
select halo pairs each of which consists of one halo from the 
low-mass bin and one halo from the high-mass bin, and then 
measure the EE correlations of these halo pairs. Figure [S] 
plots ric{r) at z = as solid dots. As can be seen, there ex- 
ist significant cross-correlations between the two mass bins. 
Note that the EE cross-correlation is in fact stronger than 
the EE auto-correlation of the low-mass halos (M < Mc) 
but weaker than the EE auto-correlation of the high-mass 
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Table 2. The EE correlations of the major axes of the low-mass 
and the high-mass halos. 
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Figure 4. Same as Fig. [3] but for the case of the halo minor axes. 
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Figure 5. The EE correlations of the halo major axes measured 
from the low-mass (M < Mc) and the high-mass (M > Mc) bins 
at 2 = in the top and bottom panel, respectively. The mass 
threshold Mc = 1416 is in unit of 10^° h-'^ Mq. 



halos (M > Afc). This result suggests that the anisotropic 
merging and infall along filaments tend to increase not only 
the EE correlations on the same mass scale but also the EE 
cross-correlations between different mass scales. 



3.4 Fitting formula 

According to the first order linear model, the EE correlations 
are induced by the spatial correlations of the primordial tidal 



log[r/(/i-lMpc)] 


rii{r) X 10^ 


r)i(r) X 10^ 




low- mass 


high-mass 


0.15 


0.78 ± 0.95 


6.28 ±2.98 


0.45 


0.74 ±0.44 


1.93 ±0.79 


0.75 


0.38 ± 0.25 


1.80 ±0.68 


1.05 


0.14 ±0.15 


0.74 ± 0.42 


1.35 


0.05 ± 0.05 


0.27 ±0.26 


1.65 


0.01 ± 0.03 


0.09 ±0.07 



10-2- 



10- 




r [Mpc/h] 

Figure 6. The EE cross-correlations of the major axes between 
the high and the low mass halos at z = 0. 



field and can be approximated in the linear regime by a 
quadr atic scaling of the linear correlation function as r)(r) oc 
^^(r) (jCatelan et al.ll200ll : lLee"fc Penll200ll : iHirata fc SeliakI 
[2OOJ). Since ^^(r) decreases rapidly with separation distance 
r, the linear model predicts that the EE correlations exist 
only between close pairs with r less than a few Mpc. In other 
words, the large scale EE correlations cannot be described 
by t he linear mod e l. 

iHui fc Zhanid (|2002l ) pointed out that the growth of the 
non-Gaussianity in the density field should cause the EE 
correlations to scale linearly with ^(r). It indicates that the 
non-Gaussianity tends to increase the la rge-scale EE corre- 
lations. Although iHirata fc Selial3 (|2004l ) discussed that the 
assumption of ri{r) oc ^(r) is valid only in the linear regime, 
we use here the following fitting formula for rj{r): 



■n{r) ~ a-Ui-r)- 



(3) 



Here a is a fitting parameter, representing the amplitude of 
the EE correlation, whose value is to be determined empiri- 
cally. Since ri{r) is always positive and does not exceed 2/3, 
the correlation parameter a is expected to be in the range 
of [0, 2/3]. 
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Table 3. Redshift (2), halo mean mass (M) in unit of 
10^^h~^MQ, number of halos (Af^), and the best-fit values of 
a for the EE correlations of the halo major axes. 



z 


M 

[IQIO/i-IMq] 


Nh 


a X 10^ 





1638.15 


121773 


1.12 ±0.06 


0.5 


1118.82 


131505 


2.45 ±0.11 


1 


810.92 


125363 


2.64 ±0.23 


2 


448.33 


73514 


3.01 ±0.39 



In eq. ((Sjl, £,(r)A is the rescaled two-point correlation 
function of the linear density field, defined as 



f P(fc)[(sin kr)/kr]W^{k;M) d^k 
/ P(fc)W2(fc;M) d^k ' 



(4) 



which satisfies Ca(0) = 1. Here, P{k) is the linear power 
spectrum, and W{k; M) is the top-hat spherical filter corre- 
sponding to the mass scale M. We have employed the ap- 
proximate formula given by iBardeen et all (Il986l ) for the 
ACDM power spectrum, using the same values of the cosmo- 
logical parameters that were used for the Millennium Run 
simulation. For the shape factor F of the power spectrum 
parameterization, we adopted F = ilmh. For the smoothing 
mass scale M in eq. Q, we use the mean mass averaged 
over the selected FOF halos. 

We fit the numerical results obtained in Section 3.2 
to eq. ((3)1 by adjusting the para meter, a with the help of 
the maximum likelihood method l|Barlowlll99ll ). Basically, 
it amounts to finding the minimum of the function de- 
fined as 



(5) 



where rji is the numerical data point at the i-th radial bin, 
ri, r]{ri;a) represents the fitting model (eq. |3} calculated at 
Ti, and {C~j^) is the inverse of the covariance matrix, iCij), 
whose component is calculated a s the ensemble average over 
the 8 realizations l)Barlowlll99lh : 



Ci 



(iVi - VOi){Vj - V03)), 



(6) 



where 7701 represents the mean rji obtained from the whole 
simulation box. It is worth mentioning here that the func- 
tion is expressed in terms of the inverse covariance matrix, 
(Cij), given that there exist non- negligible correlations be- 
tween radial bins at distance scales larger than 5/i~^Mpc. 
In our case the number of realizations is larger than that of 
radial bins, the x^ function defined in terms of the inverse 
cova riance matrix should be useful to find the best-fit value 
of a (|Hartlap et al.ll2007f) . 

The uncertainty in the measurement of a is calcu- 
lated as the curvature of the x^ function at the minimum 
l|Bevington fc Robinsonll 19961 ): 



a2 2 

d X 



(7) 



The fitting results are summarized in Table|3]which lists the 
mean mass M in unit of lQ^'^h~^MQ, the number of halos 
Nh, and the best-fit values of a at four different redshifts. 
Note that the value of a deviates from zero at all redshifts 
and decreases monotonically as z decreases. These results 



Table 4. Mass bin, halo mean mass (M) in unit of 10^^h~'^MQ, 
and the best-fit values of a for the EE correlations of the halo 
major axes. 



bin 


M 


a X 10^ 




[lOi°/i-iM0] 




low-mass 


545.7 


0.86 ±0.32 


high-mass 


4915.6 


2.75 ±0.39 



quantify well how strong the EE correlations are and how 
they evolve with redshifts. The fitting models with these 
best-fit-values of a are plotted as solid line in Fig. [S] As can 
be seen, the fitting models are in good agreement with the 
numerical results (solid dots) at all redshifts. For compari- 
son, we also try to fit the numerical results of the EE corre- 
lations of the halo major axes at z = to the linear model 
proportional to ^^(r). The linear model with the best-fit am- 
plitude is shown as dashed line in the bottom-right panel of 
FigO As can be seen, the linear model drops with r too 
rapidly to fit the large-scale EE correlations. 

We also fit the EE correlations of the halo minor axes 
measured at four redshifts and the EE correlations of the 
major axes measured from two different mass bins at 2: = 
to eq. @ and plot the results as solid lines in Figs. |3] and 
[5l The EE cross-correlation of the halo major axes between 
different mass bins at 2: = is similarly modeled as: 

nc{r)^a,^c{r). (8) 

Here Oc is a fitting parameter and (,c{r) is defined as 

/ P{k)[smkr/kr]W(k;lVh)W{k;M2)d^k 



J P{k)W{k;Mi)W{k;M2)d»k 



(9) 



where Mi and M2 represent the mean mass averaged over 
the low-mass and the high-mass bin, respectively. The best- 
fit value of Uc is determined similarly by minimizing the 
function. The fitting result for the EE cross-correlation is 
plotted as solid line in Fig. [S] As can be seen, the agreements 
between the numerical results and the fitting models are 
quite good for all cases. 



4 THE ELLIPTICITY-DIRECTION CROSS 
CORRELATIONS OF HALOS 

4.1 Definition 

As mentioned in Section 1, another important correlation 
function for quantifying the cosmic web is the ED cross- 
correlation between the halo ellipticities and the large-scale 
density field. If the halo ellipticities are induced by the 
anisotropic infall and merging along the local filaments, then 
the orientations of the halo major axes must be preferentially 
aligned with the directions to the neighboring halos. This ef- 
fect can be measured in terms of the ED cross-correlations 
between the halo orientations and the location of halo neigh- 
bors. We employ the following definition of the ED cross- 
correlations as 

..(r) = (|e(x).f(x)|^)-i, (10) 

where r = r/r is a unit vector in the direction to a neigh- 
boring halo at separation distance of r. 
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Figure 7. The ED correlations of the halo major axes (solid 
dots) at 2 = 0, 0.5, 1 and 2. The errors represent the standard 
deviation between realizations but do not include non-negligible 
correlations between radial bins. In each panel the solid line rep- 
resents the best-fit model with two parameters (eq. 111^ . In the 
bottom-right panel the dotted and dashed line correspond to the 
best-fit results based on the model with one parameter propor- 
tional to ^ and S,^^^, respectively. 

From the halo catalogs of the Millennium Run simula- 
tion at 2 = 0, 0.5, 1 and 2, we have measured the ED cross- 
correlations of the halo major (a;i), intermediate (wn) and 
minor axes (win) as a function of the comoving distance 
r between the halo centers. Basically, for each halo in the 
Millennium data at a given redshift, we find the direction 
to its neighbor halo and measure the separation distance r, 
and calculate the squares of the dot products of the nor- 
malized eigenvector with the unit vector in the direction to 
the neighbor halo. And then, we bin the radial distance r 
and calculate the mean values of |e(x) • f (x)p averaged over 
those halos whose distances to the neighbor halos belong to 
a given bin, subtracting 1/3 from it. We perform this pro- 
cedure at z = 2, 1, 0.5 and 0. 

4.2 Evolution with redshift 

Figure [7] plots lji at redshifts z = 0, 0.5, 1 and 2 in the 
top-right, top-left, bottom-right and bottom-left panel, re- 
spectively. The errors represent the standard deviation be- 
tween 8 realizations. Since we are mainly interested in the 
cross-correlations between the halo principal axes and the 
large scale density field, we focus on separation scales greater 
than l/i~^Mpc. As can be seen, the ED cross-correlations of 
the halo major axes also decrease as z decreases at all dis- 
tance scales. Note also that the ED cross-correlations are 
much stronger than the EE correlations shown in Fig. [S] 
The ED signal is statistically significant even at distances 
out to 50 /i"^ Mpc. 

The ED cross-correlations of the intermediate and mi- 
nor axes of halos, uju and ujiu, are plotted in Figs. [8]and|9l 



Figure 8. Same as Fig.[7]but for the case of the halo intermediate 
axes. 
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Figure 9. Same as Fig. [7]but for the case of the halo minor axes. 



respectively. As expected, the intermediate and minor axes 
are anti- correlated with the directions to neighboring halos, 
and the degree of the anti-correlation is stronger for the mi- 
nor axes. In fact, the ED anti-correlations of the halo minor 
axes are almost as strong as the ED correlations of the halo 
major axes. These results demonstrate clearly that the halo 
major axes preferentially point in the directions where the 
local density stays high, hence this gives a quantitative mea- 
sure for the filamentary distribution of the halos in the cos- 
mic web. The numerical results on uii, lju and uim measured 
at z=0 are summarized in Table [5] 
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Table 5. Numerical results for the ED cross-correlations of the 
halo major (tiJi), intermediate (i^ii), and minor (ljiii) axes in log- 
arithmic scale measured at z = 0. 



log[r/(fe-iMpc)] 


uji{r) X 10^ 


-u)li(r) X 10^ 


—u)iii{r) X 10^ 


0.375 


6.85 ±0.35 


1.19 ± 0.42 


5.67 ± 0.40 


0.625 


5.37 ±0.26 


1.04 ±0.16 


4.33 ± 0.22 


0.875 


3.42 ±0.17 


0.73 ± 0.12 


2.68 ± 0.12 


1.125 


1.80 ±0.11 


0.36 ± 0.08 


1.43 ± 0.07 


1.375 


0.88 ± 0.05 


0.20 ± 0.03 


0.68 ± 0.06 


1.625 


0.37 ±0.04 


0.11 ± 0.03 


0.26 ± 0.04 
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Figure 10. The ED cross-correlations of the halo major axes 
measured from the low-mass (M < Mc) and the high-mass (M > 
Mc) bins at z = in the top and the bottom panel, respectively. 
The mass threshold Mc = 1416 is in unit of 10^° H-'^Mq. . 

4.3 Variation with mass 

To explore how the ED correlation changes with halo mass, 
we measure the ED correlation of the halo major axes at two 
different mass bins with the mass threshold Afc = 1-42 x 
10"/i"^Mq at z = 0. . When finding the neighbors, we 
consider all halos, no matter what mass the neighbor halos 
have. Figure [TO] plots the ED correlations of the halo major 
axes at two different mass bins at 2 = 0. As can be seen, 
the ED correlations increase as the halo mass increases, just 
like the EE correlations, which suggests that the anisotropic 
merging contributes significantly to the ED correlations. 

4.4 Fitting formula 

Given the observed fact that the EE correlations scale lin- 
early with the linear density two-point correlation function 
(Section 3.4), we have also tried in vain to model the ED cor- 
relations as a linear scaling of the density correlation func- 
tion. But, it has turned out that this simple model fails in 
providing good fits to the large-scale ED cross correlations 
(see Fig. [Tjl. To improve the fitting result, we need a model 



which decreases with r more slowly than S,{r). We suggest 
the following empirical formula: 

i^{r)^hiA{r)+b2ef{r) (11) 

where the two parameters 61 and 62 lie in the range of 
[—1/3, 2/3]. The second term proportional to ^^/^ is in- 
cluded to fit the large-scale ED correlations. For the halo 
major-axis, both of bi and 62 will have positive values, while 
for the halo minor and intermediate axes, they will be neg- 
ative. We also expect that bi and 62 will have larger values 
than a, since the ED cross-correlation is a more direct mea- 
sure of the filamentary distribution of dark matter halos. 

We fit the numerical results on the ED correlations ob- 
tained in Sections 4.2-4.3 to eq. and determine the best- 
fit values of 61 and &2 by minimizing the function: 

= ['^^ ~ i^{r^; hi, b2)]Cr/[ujj - (.^(rj; 61, 62)], (12) 

where uii is the numerical data point at the z-th distance bin, 
ri, and aj(ri; fei, 62) is the fitting model at r^. The covariance 
matrix, (Cij), is calculated as 

Cij = {{uJi - LUOi){u>j - l^Oj)), (13) 

where ojoi is the mean u>i obtained from the whole simulation 
box. To calculate errors in the measurement of 61 and 62, 
we first construct a 2 x 2 curvature matrix defined as 

with i,j = {1,2}. The errors are calculated as t he diago- 
nal c omponents of the inverse curvature matrix (jDodelsonl 
l2003h : 

4 = (^"')22- (15) 

Table [5] lists the best-fit values of 61 and 62 for the ED 
correlations of the halo major axes measured at 2: = 0, 0.5, 1 
and 2. The best-fit value of the parameter 62 is larger than 
that of hi at every redshift, indicating that the second term 
in eq. (Illf) dominates. The values of the two parameters 
decrease monotonically as z decreases, just like the EE cor- 
relations. Our results provide a quantitative description of 
the evolution of the ED correlations and its scaling with 
distance. 

The fitting results on the ED correlations of lji, uju and 
ojiii at z = with the best-fit values of 61 and &2 are plot- 
ted as solid lines in Figs. [T] [SI and [§] respectively. As can 
be seen, the agreements between the fitting models and the 
numerical results are quite good at all redshifts. For compar- 
ison, we also fit the numerical results of the ED correlations 
of the halo major axes at 2 = to two different models 
proportional to ^ and which are plotted in the bottom 
right panel as dotted and dashed line, respectively. As can 
be seen, the model proportional to ^ decreases with r too 
rapidly to fit the numerical results. Meanwhile the ^^''^ de- 
creases slowly with r but it alone still does not agree with 
the numerical results as well as eq. (Illf) . 

Figure[lO]plots the fitting results on the ED correlations 
of the major axes of the low-mass and high-mass halos at z — 
in the top and bottom panel, respectively. Table[7]lists the 
best-fit values of fei and 62 for the two cases. The values of 61 
and 62 are higher for the high-mass halos than for the low- 
mass halos. Note that although the fitting models work fairly 
well, they seem to deviate from the numerical results by 
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Table 6. The best-fit values of the two correlation parameters 
for the ED correlations of the halo major axes at four different 
redshifts. 



z 


bi X 10^ 


62 X 10^ 





3.06 ±0.20 


4.61 ±0.11 


0.5 


4.54 ±0.21 


5.55 ±0.12 


1 


5.73 ± 0.23 


7.25 ±0.10 


2 


7.57 ±0.32 


11.64 ±0.16 



Table 7. The best-fit parameters of the ED cross-correlations 
measured from two different mass bins at 2 = 0. 



mass bin 


M 


bi X 10^ 


62 X 10^ 




[IOIO/1-IM0] 






low- mass 


545.7 


1.39 ±0.36 


4.02 ± 0.14 


high-mass 


4914.35 


4.94 ±0.29 


4.39 ± 0.16 



more than Sa^ around 30 Mpc. We think that it reflects 
the failure of the assumption that the ED correlations can be 
described in terms of the linear density correlation function. 



5 SUMMARY AND DISCUSSION 

In this work, by analyzing the halo data and the semi- 
analytic galaxy catalog from the Millennium simulations 
at z — 0, 1, 0.5 and 2, we have measured the ellipticity- 
ellipticity (EE) correlations. The correlations are close to 
0.01 at a distance of l/i~^Mpc and remain significant at 
distances out to 10 Mpc. The EE correlations are found 
to be strongest for the case of the halo major axes, and 
weakest for the case of the intermediate axes. 

We have found that the EE correlations of all three axes 
decrease as z decreases. This might be due to the growth 
of secondary filaments at low redshifts and the beginning 
'freeze-out' of structure growth in ACDM, which plays a role 
in randomizing the halo ellipticities. It has been also found 
that the EE correlation function exhibits a strong depen- 
dence on halo mass. It increases as the mass scale increases, 
which might be due to the dominant filamentary merging 
of halos on large scales. We have also calculated EE cross- 
correlations between halos belonging to different mass bins. 
Our results have shown that the EE cross-correlations be- 
tween neighboring mass bins exist at a statistically signifi- 
cant level as well. 

We have modeled the EE correlation function as a linear 
scaling of the linear density two-point correlation function 
5(r), which is characterized by one fitting parameter, a. The 
value of a represents the amplitude of the EE correlations, 
quantifying its scaling with mass and redshift. The fitting 
model with the best-fit value of a has been shown to agree 
with the numerical results quite well at all redshifts and all 
mass bins. 

We have also measured the cross-correlations between 
the halo principal axes and the directions to neighboring 
halos (ED) by using the same numerical data, and found 
that the ED cross-correlations are much stronger than the 
EE correlations, at all distances. Remarkably, they are de- 



tected even at distances out to 50/i~^Mpc at a statistically 
significant level. Just like the EE correlations, the ED cross- 
correlations are found to decrease as z decreases and increase 
as the halo mass M increases, suggesting a dominant role 
of anisotropic merging and infall of matter in establishing 
these correlations. The intermediate and the minor axes of 
the halos have turned out to be anti-correlated with the di- 
rections to the neighboring halos, which is consistent with 
alignments of the halos shapes with the orientations of the 
local filament. 

The ED cross-correlations are, however, found to be 
poorly fitted by a linear scaling of £,{r). The ED cross- 
correlations decrease with distance much less rapidly than 
^(r). To account for the slow decrease of ED cross- 
correlations with distance, we include an additional term 
proportional to C^''^(r) in the fitting model which is then 
characterized by two fitting parameters bi and &2. Thetwo 
parameters represent the amplitudes of the two terms of 
the ED cross-correlations proportional to ^(r) and ^^^^{r), 
respectively. This fitting formula has been shown to agree 
with the numerical results quite well at all redshifts and at 
all mass bins. 

Nonetheless, it is worth mentioning here that our fitting 
formula for the ED cross-correlations is not a physical model. 
It is purely empirical, obtained through comparison with the 
numerical results. It has yet to be understood why the ED 
cross-correlations scale as described by our fitting formula. 
At any rate, we believe that our fitting formula may be useful 
in the future when the ED cross-correlations can be modeled 
by a fundamental theory. 

The EE correlations and the ED cross-correlations that 
we have measured here provide a useful tool to statistically 
characterize the anisotropy and the relevant scales of the 
cosmic web. It will be interesting to compare the results we 
obtained here for the ACDM cosmology with observational 
data from large galaxy redshift surveys. The comparison of 
our numerical results with observational data, however, will 
require a modelling of the redshift space distortions as well 
as of two dimensional projection effects, given that in real 
observations what can be usually measured is the two di- 
mensional projected major axes of the galaxies in redshift 
space. In future work, we plan to model these two effects 
on the EE and ED correlations and compare the numerical 
results with real observational data. 

Another important application of our results lies in 
studies of weak gravitational lensing. The issue of a po- 
tential cross-correlation between galaxy ellipticities and the 
weak gravitati onal lensing shear (GI c ross-correlations) was 
first raised by iHirata &: SeliakI (|2004l '). They claimed that 
if such GI cross-correlations exist, then they would af- 
fect the weak lensing signal as another systematic con- 
taminant whose effect is hard to control. The GI cross- 
correlations are expected to occur primarily due to the ED 
cross-correlations: If the intrinsic ellipticities of the galaxies 
are cross-correlated with the surrounding large-scale den- 
sity field, then it will in turn lead to a cross-correlation be- 
tween the gravitational lensing shear and the galaxy ellip- 
ticities. Recent observations indeed have reported detections 
of the GI correlation signals in low-redshift g alaxy surveys 
(|Mandelbaum et al.l [20061 : iHirata et al"]|2007l ). To assess a 
possible systematic contamination of weak lensing due to the 
GI cross-correlations, it will be important first to examine 



10 J. Lee, V. Springel, U.L. Pen and G. Lemson 



the relation between the observed GI cross-correlations and 
the ED cross-correlations of the cosmic web. This work will 
require incorporating a model for how the galaxy shapes are 
aligned relative to the dark matter tlfevmans et al.i 
Our future work is in this direction. 
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